POSITIVITY AND LOWER BOUNDS TO THE DECAY 
OF THE ATOMIC ONE-ELECTRON DENSITY 
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Abstract. We investigate properties of the spherically averaged 
atomic one-electron density p(r) . For a p which stems from a phys- 
ical ground state we prove that p > 0. We also give exponentially 
decreasing lower bounds to p in the case when the eigenvalue is 



qq ' below the corresponding essential spectrum. 

> 

^ ■ 1. Introduction and results 
O ■ 

00 ■ Let H be the non-relativistic Schrodinger operator of an iV-electron 

atom with nuclear charge Z in the fixed nucleus approximation, 

o : N f Z \ 1 

Mh' j=l 1 n l<i<j<N 1 

^ ! Here the x,- = (xj i, a; 3 - 2, 2^,3) G K 3 , j = 1, • • • , N, denote the positions 



of the electrons, and the Aj are the associated Laplacians so that A = 
J2f=i Aj is the 3iV-dimensional Laplacian. Let x = (xi, x%, . . . , xn) G 
; ]R 3Ar and let V = (Vi,...,Vjv) denote the 3iV-dimensional gradient 

^ ■ operator. We write 

S ; = —A + V (1.2) 

where V is the multiplicative potential 

V Z 



v ^ = J2~n + E 1 1- ( L3 ) 



j=l 1 - 71 l<i<j<N 1 Jl 



d37V> 



The operator if is selfadjoint with operator domain T>(H) = W 2,2 ( 
and form domain Q(H) = W 1 ' 2 ^) 0. 

For an eigenfunction ip G L 2 (M. 3N ) of H, with eigenvalue J5, i.e., 

Hip = Ei/j, (1.4) 
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we associate the one-electron density p G /^(R 3 ). It is defined by 

N N „ 

Pip) = J^Pjix) = J2 \ij)(x,Xj)\ 2 dxj , (1.5) 

where we use the notation 

Xj := (xi, . . . , xj-!, x j+ i, . . . , x N ) G R 3Af ~ 3 , 

and 

dStj := dx\ . . . dxj-idxj + \ . . . dxjy, 

and, by abuse of notation, we identify (xi, . . . ,Xj_i,x,Xj + i, . . . ,xn) 
with (x, Xj). The spherically averaged density p G L 1 (R + ; r 2 dr) is then 
defined by 

N N 

P(r) = Y] Pjif) = J2 / Pj( ruJ ) duJ , C 1 - 6 ) 

where r = u = x/\x\ G S> 2 for x = ruj G M 3 . 

We assume throughout when studying p that E and ^ in f)1.4j) are 
such that there exist constants C ,7 > such that 

|^(x)| < C e~ 7|x| for all x G R 3N . (1.7) 

The a priori estimate [TJ Theorem 1.2] (see also [7[ Remark 1.7]) and 
(jl.7|) imply the existence of constants Ci,7i > such that 

|V^(x)| < Ci e~ 7l|x| for almost all x G R 3N . (1.8) 

Remark 1.1. Since ■?/> is continuous (see [S]), (jl.7|) is only an assump- 
tion on the behaviour at infinity. For references on the exponential 
decay of eigenfunctions, see e.g. [3] and [TJ]. The proofs of our re- 
sults rely (if not indicated otherwise) on some kind of decay-rate for 
■0; exponential decay is not essential, but assumed for convenience. 
Note that (|1.7|) and (jl.8J) imply that p is Lipschitz continuous in R 3 
by Lebesgue's theorem on dominated convergence. This on the other 
hand implies that p is Lipschitz continuous in [0, oo). 

Since we are interested in atoms (in this non-relativistic description 
with fixed nucleus) we have to take into account that electrons are 
fermions. We shall work in the spin-independent description. That is, 
we split N such that 

N = N 1 + N 2 , N u N 2 >0, 
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and proceed as follows: We associate Sobolev-spaces to this splititng. 
Let iS(R 3iV ) be the space of Schwartz-functions, and define 

S NuN2 (R 3N ) = {ip G <S(R 37V ) | (p(x 1 ,x 2 ,x 2 , . . .,x Nl ,x Nl+1 , . . .,x N ) 

is antisymmetric with respect to the first Ni coordinates 

and antisymmetric in the remaining N2 coordinates. } 

Therefore, for instance, 

j X% , . . . , Xj , . . . , X jVj , . . . , Xn ) 
= — (f{Xi, . . . , Xj, . . . , Xi, . . . , Xn-l , ■ ■ ■ , Xjy). 

Similarily (p changes sign if we interchange the coordinates of two elec- 
trons which belong to the other group of N2 electrons which are labeled 
with i = Ni + 1, . . . , N. Note that in physical terms this require- 
ment means that the total spin is ±|iV2 — N±\/2. Define finally the 
Sobolev-spaces N2 (M. 3N ) as the closure in the M^ p,l3 (]R 3Ar )-norm of 

Let H Ni ^n 2 be the atomic iV-electron Schrodinger operator defined 
by (jl.lj) . restricted to functions with the above symmetry. Then Hn 1 ,n 2 
has operator domain V(Hn 1 ,n 2 ) — Wivf n 2 (^- 3N ) an d f° rm domain 
Q{H NltN2 ) = W^^R 3 ^). We denote E NljN2 the infimum of its spec- 
trum (when this is an eigenvalue), and call it the ground state energy. 
A corresponding eigenfunction ip = ipN 1 ,N 2 is called a ground state. E 
will henceforth denote any eigenvalue. 

Here is the first of our main results. 

Theorem 1.2. Let ip be a ground state of H Nl N2 , i.e., H NltN2 ip = 
En!,n 2 iP, and let p be the associated spherically averaged density defined 

by (jlSD-jEEIl. 

Then 

p{r) > for allr G [0, 00) . (1.9) 

Remark 1.3. 

(i) At the origin we derive an explicit, positive lower bound to the 
density (see ([2~T3]^ 

2 p4 .. N 

''^kw • withpH|y>,4 (L10) 

(ii) Note that the choice of anti- symmetric in both groups of coor- 
dinates in the definition of iSa^a^R 3 ^) is, in fact, not essential. 
One could consider functions symmetric in each group of co- 
ordinates. In fact, the theorem holds for any combination of 
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symmetric/anti-symmetric. This will be clear from the proof. 
In particular, with Ni = N and symmetric (N 2 = 0), one gets 
the known result for the absolute (bosonic) ground state. 

(iii) We would expect that the non-averaged density p of a ground 
state of H NljN2 does not vanish either. Also, the one-electron 
density p associated to fermionic ground states of molecules 
should be strictly positive. However, these are much harder 
problems and they will not be addressed here. 

(iv) It is not clear what to expect for excited states. Consider for 
instance a two-electron atom with no interelectronic repulsion; 
that is, H = —Ax — A 2 — Z/\xi\ — Z/|x 2 |. Let 4>u i — 1, 2, 
be linearly independent eigenf unctions of the three-dimensional 
one-electron operator —A — Z/\x\ satisfying </>i(0) = 02(0) = 0. 
Then i])(x\,X2) = <pi(xi)<p2{x2) — <fii (^2)^2(^1) is an eigenfunc- 
tion of if2,o such that the associated density satisfies p(0) = 0. 
However, it is not clear whether or not p still vanishes once the 
interelectronic repulsion is turned on. 

Our next result on the density is in the spirit of [3]. 

Theorem 1.4. Let ip be an eigenfunction of Hn 1 ,n 2 with eigenvalue 
E and let p be the associated spherically averaged density defined by 
f01)- ffH)J) . Define 

a = sup{a|e Q|x| ^GL 2 (M 3JV )}. (1.11) 



Then 



luii sup ) < -2of . (1-12) 



R- 



If furthermore E < inf o"ess(-f^7Vi,Af 2 ); then also 

l imin f (^M) >-2VAW (1-13) 

Remark 1.5. One can make these bounds more explicit using ^j; in 
fact, 

al<\E\. (1.14) 

To see this, we use Theorems 1.1 and 1.2 in 0]. The set of thresholds 
T(H) is defined as the closure of the set of eigenvalues of subsystems, 
i.e., the corresponding ionized systems. We have, according to jl], 

al + E e T(H) and T(H) C (-00, 0] , (1.15) 

so that 

al < sup T(H) -E= \E\. (1.16) 
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We shall discuss in Section Kill why the difference beween the upper 
and lower bounds (|1.12j) . (|1.13|) is to be expected. 

Remark 1.6. It will be clear from the proof of Theorem 11.41 (see 
Section EJ) that the result holds for more general iV-body operators. 
That is, we could replace V in fjl.3|) by any multiplication operator of 
the form ^2 u v u (Il u x.), x G M. 3N , where U u is the orthogonal projec- 
tion of M. 3N on a <i„-dimensional subspace, and the v v are real-valued 
functions on M. d " satisfying that v u (—Ad„ + 1) _1 and (— + l)^ 1 (y • 
Vv u (y))(— A dv + 1) -1 are compact as operators on L 2 (M. d "). Here —A dv 
is the usual Laplace operator on M. " , and y G IR " (compare with (1.2), 
(1.3) and (1.4) in jlj). 

The above theorem gives upper and lower exponential bounds on p 
near infinity. Combined with Theorem 11.21 this implies (by continuity 
of p, see Remark 1 1.1 Jl global lower exponential bounds in the case of a 
ground state. We state this explicitely in the next corollary. 

Corollary 1.7. Let be an eigenf unction of H^ 1) n 2 with eigenvalue E, 
let cto be as in (jl.llj) . and let p be the associated spherically averaged 
density defined by ()1.5|) - (jl.6j) . 

If E < inf c ess (H NljN2 ) , then for all a > a there exists r > and 
c = c(a, r ) > such that 

p{r) > ce- 2VWar for all r > r . (1.17) 

If furthermore E = En u n 2 (the ground state energy), then p,17j) 
holds with r = 0. 

More detailed results than ()1.12J) on upper bounds to p are known. 
For completeness, we include the following classical result; see 0. 

Proposition 1.8. Letip be an eigenfunction ofH^ 1 ^ 2 with eigenvalue 
E, satisfying 

e:=mfa ess {H NuN2 )-E>0, (1.18) 

and let p be the associated density defined by (J1.5J) . 

Then for all r > there exists a constant C = C(r ) > such that 

p(x) < C \x\ ^ e-V^I f or a u x e with \x\ > r . (1.19) 
Remark 1.9. 

(i) The infimum of the essential spectrum of ifjVi,7v 2 is character- 
ized by the HVZ-theorem |111 Theorem XII. 17], which takes 
symmetry into account, in particular the fact that we consider 
fermions. Hence, 

e > min [E Nl _ 1)N2 , .EVi^-i} ~ E. 
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(ii) Note that the bound (J1.19j) can be asymptotically sharp as has 
been shown for the ground state density of the Helium-like atom 

where the physical ground state eigenfunction can be chosen 
positive. 

(iii) If we consider a bosonic ground state, Hip = E Q ip, for an atomic 
Hamiltonian which happens to have Eq = inf a ess (H) then the 
associated density can decay like e _/3 ^ for some suitable (3 up 
to some sub-exponential factors; see jHj. 

2. Proof of positivity 

Proof of Theorem \l. 6 A Assume first for contradiction that 

p{ r o) — for some r > 0. (2.1) 

This implies that Pj(ro) = for all j = 1, . . . , N, and therefore, still 
for all j = 1, . . . , N, that 



x. 



12 - if x G Nj{r ) = {x G R 3N \ \x,\ = r }. 



Here we used fll.5|) . ()1.6|) . and the continuity of if). This means that 

V>(x) = for x G N(r Q ), 

where 

N N 

iV(r ) = |J N 3 (r ) = \J {x G R 3N \ \ Xj \ = r }. 

j=l 3=1 

Let 

O = ri (r ) = {x G R m I max \xj\ < r \. 

j=l,...,N 

Then Oo is an open bounded subset of R 3N satisfying 

x G fl => "Px G Oo 

for every permutation V G & N of the electron coordinates (that is, 
of the iV-tuple (1, . . . , N)). This means that the following space is 
non-trivial (7^ {0}): 

Q Nl , N2 (n ) = {fe <' 2 (O ) j 3F G Q(H NliN2 ) such that F\ no = /}. 

In fact, by the above, ip = on dflo. Therefore, the restriction ipn : = 
ip\a of ip G Q{H NlyN2 ) to O belongs to Qjvi,jv 2 (Oq), and clearly 

(Vno,^n > = ^,^11^0 f- (2-2) 



POSITIVITY AND DECAY OF THE ELECTRON DENSITY 



7 



We now claim that we have the strict inequality 

<peQN lt N 2 {n ) \\(p\\ 2 ■ <peQ(H NltN2 ) \\(f\\ 2 

(2-3) 

Indeed, assume that we have equality in (|2.3|) . Then, by the variational 
characterization of the ground state (also in a symmetry subspace), the 
eigenfunction which minimizes the LHS of ()2.3|) . and which we extend 
to M. 3N by setting it identically equal to zero outside Qq, will have to 
be an eigenfunction in all of M 3Ar also. This is a contradiction to unique 
continuation (see [TT] Theorem XIII. 57]). 

Now, (|2.2|) and (|2.3|) imply that i()\n = 0. By unique continuation, 
■0 = (since r > 0). This is a contradiction, and therefore settles the 
case when r > (see (|2.1j) ). 

We still have to show that 

p(0) = 47rp(0) > 0. (2.4) 

Here we explicitely use the Coulombic nature of the potential V (see 
(|1.3j) ). Recall that x = (x±, . . . , xn) G R 3JV , with Xj = (xj t x, Xj^, ^,3) G 
M 3 . The gradient with respect to Xj is denoted V,-. 
Let, for «6l, 

N 

F [a) = + ax ^ = i F i a) ,Ft\ F 3 (Q) ) (2.5) 

3=1 

with F { k a) = J2f=i + oix jtk ip) , k = 1, 2, 3. Using that ^ G £>(F) = 
VY 2 ' 2 (]R 3Ar ), and the exponential decay flTTD and ljTBjl . we get that 
G H /1,2 (lR 3Ar ). Then the following variational expression is well- 
defined: 

R («) = E <^ - E Nl , N2 )F^). (2.6) 



fc=i 



Note that the F^ are obtained by applying the operator 

N q 
, dx 3,k 

to ip, which does not change symmetry properties with respect to per- 
mutations of the electron coordinates. Hence has the same per- 

mutational properties as ip itself, and so G Q(Hn u n 2 ), k = 1,2, 3. 
Therefore, by the variational characterization of the ground state ip 
(also in a symmetry subspace), we have 

R(a) > for all a G R. (2.7) 
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We have (since H is self-adjoint) that 

3 N 



k=l ij=l 



dip dip 



, dip 



(H - E Nl , N2 )(x i>h ip)), 



Ci,j,k = ((xj, k ip), {H - E NuN2 )(x itk ip)). (2.8) 
First note that, since Hip = E NltN2 ip, ([ ; ] denoting the commutator) 



(H-E 



dip 



N lt N 2/ 



[H-E 







dx i:k 

and so, by partial integration, 

3 N 3 AT 



N u N 2 j 



dx 



i,k 



k=l i,j=l k=l i,j=l J ' 



dV 



3 N 



k=l i,j=l 



k d%i,k 

d 2 V 

dx jyk dx i>h 



)*) 



Strictly speaking, V and ip are not smooth enough for this and the 
following computation to be but formal. However, regularizing V and 
using the exponential decay (jl.7|) - (jl.8j) . and the continuity of p (see 
Remark II. 1)1 . one easily justifies this. 
Note that (for V, see (Oil ) 

t 9 d v. /, . n 

(e^ + a^ )ilx, - Xil ) = °' 



E 

fc=l 



dxj^dxi^ 



where Sij is Kronecker's delta, and S is the delta - function. This way, 



3 N N 
k=l i,j=l j=l 



2nZ 



(x)| 2 5(|^|)dx 



N 



27rZ^p j (0) = 2ttZ P (0). 



(2.9) 



POSITIVITY AND DECAY OF THE ELECTRON DENSITY 9 

Secondly, again since Hip = E^^if), 

(H - E NljN2 )(x itk iP) = [(H - E NuN2 ); x hk }^ = - 2^-. (2.10) 



Therefore 

dip \ ^ dip 



Y Y hi ^ k - ~ 2 Y(Y a^T' Y -fa—) 

k=l i,]=l k=l 3=1 J ' 1=1 

II N 2 

= ~ 2 \\Y V ^ ='-- 2p2 - ( 2 - n ) 



3=1 

Finally, using ()2.10|) and partial integration, 

3 N 3 N dip 

Y Y Ci >j> k = ~ 2 YY ( x ^> -fa-) 

k=l i,j=l k=l i,j=l 

3 N 



YY I x ^^-{^ 2 )d-> 



k=l i,j=l ' 
3 Af 

= YY 6 ^r\m 2 = mm 2 - (2.12) 

k=l i,j=l 

Combining flZSJl, (EZH), and ([2~T2jl with ((231), we find tliat 

R(a) = 2nZp(0) - 4aP 2 + 3a 2 iV||^|| 2 > for all aeK. 
Optimizing this expression in a we obtain 

j=l 

Suppose that P = 0. After Fourier transformation, this implies that 

5>i1?(p) = °, (2.14) 

3=1 

in L 2 (R 3Ar ). The equation (|2.14j) clearly implies that ip = (since it has 
support on the set { p £ R 3Ar | ^2j=iPj = ® } which has zero measure). 
Since if) is an eigenfunction, ip ^ 0, so this is a contradiction. We 
conclude that P^O and (l2~T3|) thus implies (l2~3J) . 

This settles the case r = 0, and therefore finishes the proof of The- 
orem 11.21 □ 
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3. Discussion and proof of decay 

3.1. Discussion and examples. Before proving Theorem 11.41 we 
show with some examples why the difference on the right hand sides of 
()1.12j) and ()1.13|) is natural. Define 5?o as 

5 = sup {a | e 2ar p(r) G L 1 (K+ , r 2 dr)} . (3.1) 

For the definition of a , see (jl.ll|) . 

Example 3.1 (5o = «o)- Consider first the function, 

^(xi, x 2 ) = e -*i\*i\ e -*2\^\ } Xl ,x 2 e M 3 . 

Clearly, the associated density satisfies, 

p(r) = Cl e~ 2air + c 2 e~ 2a2r , c x = vr/a 3 , c 2 = ir/af . 

Thus, 5o = min(ai, a 2 ). It is not hard to see that in this example also 
a = min(a;i,a 2 ): Clearly, a > min(ai,o; 2 ). Suppose, without loss of 
generality, that a\ < a 2 . If a > ai, then (by continuity) there exists a 
neighbourhood T (in S 5 ) of (1, 0, 0, 0, 0, 0), and e > 0, such that, for all 
(uji,uj 2 ) G r, 

a — (Xi\oji\ — a 2 \uj 2 \ > e. 

The integral of le"'"''?/'! 2 over the cone M + T C M 6 therefore clearly 
diverges. Thus, 

ao = min(ai, a 2 ). 
Example 3.2 (5 7^ « ). Consider now the function 

ijj 2 {x u x 2 ) = e - ail — l e- a21 — 1 , x u x 2 G R 3 . 

It is easy to see that ao = min(ai, a 2 ) also in this case (the definition of 
cto is invariant under an orthogonal change of coordinates). However, 
we will see that 5 = \/2 min(«i, a 2 ). We write out 

/ e 2a ^ Pl (x 1 )dx 1 = [ e^e-^^e^^dxrdx, 

JR3 J R 6 

= [ e 2a ^e- 2a ^e- 2a ^d yi dy 2 . 

Since \yi + y 2 \ < + \y 2 \, the above integral is clearly convergent 
for all a satisfying < min(a 1 ,a2). It is also easy to see (as in the 
previous example) that the integral is not convergent (on a suitable 
cone) if > min(ai, a 2 ). 
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Remark 3.3. Example 13.11 would be the correct behaviour (modulo 
polynomial forefactors) of eigenfunctions of H, if we omitted the terms 
V. Since our proof works for general iV-body opera- 
tors (see Remark ll.6|) . it is therefore clear that Example 13.21 is equally 
relevant, since Example 13.21 is obtained by using a rotation of the co- 
ordinates in Example 13.11 

The proof of Theorem 11.41 will rely on the result pfl Theorem 2.1] 
adapted for our purpose. For Ri < R 2 , denote by Q(R\, R 2 ) the spher- 
ical shell (in 



0(i?i, R 2 ) = {x € R m | Ri < |x| < R 2 }. (3.2) 

Theorem 3.4. Suppose if) is an eigenfunction of H with eigenvalue E 
and let a be defined by (I1.11J) . Suppose S(R) is a positive function 
satisfying 

hi 5(R) \ 

Then 



liminf — ^ > 0. (3.3) 



lim 4 m ( / |^(x)| 2 dx) = -2a - (3.4) 



-R^+oo R 



Q(R,R+S(R)) 



The next result, which is a special case of 01 Theorem 2.2], will 
not be used in the sequel. It is given only in order for the reader to 
be able to compare the result for the density, Theorem 11.41 with the 
corresponding result for the spherically averaged eigenfunction. 



Theorem 3.5. Suppose in addition to the hypotheses of Theorem U~J 
that E < inf a ess (H). 
Then 

lim 4 ln f / \ip(R£i)\ 2 dn) = -2oq. 

R—>+oo R V Jg3N-i J 

Before starting the proof of Theorem II A\ we state and prove a result 
similar to Theorem 13 .41 for the density p. 

Theorem 3.6. Suppose if) is an eigenfunction of H with associated 
electron density p and let ceo be defined by (jl.llj) . 
Then 

lim sup — In ( f p(r)r 2 dr) < — 2a , (3.5) 

ft^+oo R \J R / 

liminf — In ( / p(r)r 2 dr ) > -2\^Na Q . (3.6) 

i?,^+oo R \ J R J 
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Proof. The proof of ()3.5j) is a simple calculation: 

POO P 

i p(r)r 2 dr < N |^(x)| 2 rfx 

JR J {x : max j ■ |xj \>R} 

<N I \ip(x)\ 2 dx. 

J{x : |x|>R} 

Therefore, for all a < a , 

POO 

/ p(r)r 2 dr < Ne~ 2aR || e a|x| ^ || 2 2 

This clearly implies (|H.5jl . 

In order to prove ()3.6j) . we calculate 



oo 

2 



p{r)r 2 dr 

N poo p p 

= \ / r 2 dr duj \i/j(xi, . . . , XN)\ 2 S(rui — Xj) d^ 

N 

= V / |^(xi,...,Xiv)| 2 o?x 

> / \i>{xi, ■ ■ ■ ,x N )\ 2 dx. 

J {x : maxj \xj | >R} 



Using 



we get 



min \xj\ < —== < max ij , 
3 \/N ~ 3 3 



{x G K 3Af I VNR < |x|j c{xG M 3Ar I max|zi| > R\, 

j 

and so, from the above, 

POO P 

/ p(r)r 2 dr > / \tj;(xi, . . . ,x N )\ 2 dx.. 

JR J{x:VNR<\x\} 

The lower bound (|3.6|) . now follows from Theorem 13. 41 bv taking S(R) = 
+oo. □ 

3.2. Proof of Theorem 11.41 We now start the proof of Theorem II .41 
It will be based on a number of lemmas. 

We start by defining an operator Tr which maps functions from 
#1/2+6^3^ con ti nuous ly to L 2 (§ 2 x M 3iV - 3 ). T R is the restriction map 

[T R f] (u, Xl ) = f(Roo, Xl ) , R > , u G § 2 , Xl G M 3 ^ 3 . (3.7) 
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Lemma EI3 will allow us to get the upper bound ()1.12j) . 

Lemma 3.7. Suppose ip is an eigenf unction of H with eigenvalue E. 
Then there exists c > such that, for all R > 1 and all j G {1, . . . , N}, 

POO 

p j (R)<cR~ 2 pj(r)r 2 dr. (3.8) 

JR-l 

In particular, 

poo 

p(R)<cR- 2 p(r)r 2 dr. (3.9) 

J R-l 

Proof. Since (J3.9|) follows from (J3.8)) by summation over j, it clearly 
suffices to prove the latter. Without loss of generality, we will restrict 
attention to the case j '• = 1. It is well-known that the restriction map 
Tr from ()3.7|) defines a bounded map between Sobolev spaces with loss 
of a half (+e) derivative. We need some control of how the constants 
depend on R, but not the optimal regularity result, so we state and 
prove the following (elementary, not optimal) lemma. 

Lemma 3.8. The map Tr defines a bounded operator from H 1 (R 3N ) 
to L 2 (E> 2 x M 37V - 3 ) ; with the estimate 

||^i?/||L 2 (S 2 xR3JV-3) < -R _1 (||V/|| L 2( R 3iV) + ||/||l2( M 3JV)) . (3.10) 

Proof of Lemma \3. M (The proof is a repetition of the proof of 
jH Lemma 3.1]). For / e C£°(M 3Ar ), we have 

\f(Ru, xi)| 2 duj dk\ 

; J xR 3iV-3 

00 r d ... _, 2 



— — \f(ru, xi)| dudxidr 

r°° r d 

— I / — 2Re{/(ru;, X4)— f(ru, xi)} du dki dr 

POO p 

<R~ 2 / (Ifircu^^ + lVfiruj^^ydujd^dr. 

Jr j§ 2 xr 3JV - 3 

This clearly implies (|3.1U|) . from which Lemma [3.81 follows. □ 

We now finish the proof of Lemma l3~7l Let \ £ C°°(R) be monotone, 
< X < 1) sucn that 

X(t) = for t < 1/2 , x(t) = 1 for t > 1. (3.11) 

Define Yj?,( x ) — xd^il — (R ~ !))• With ijj being the eigenfunction of 
H, we have Tri/j = Tr(xriP), and therefore (using the relative form 
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boundedness of V with respect to the Laplacian Theorem X.19]) 
we get, using Lemma l3~51 

p x (R) = \\T R ip\\l2 {S 2 xM 3N- 3) < R~ 2 (xr^, (-A + l)xRtp) 
<cR- 2 ( X R^AH+l)XR^} 
= ciT 2 ^, {\{ X 2 rH + HxD + xl + I Vxi?| 2 )^) 
= cR~ 2 ^, ((E + 1)x 2 r +\V X r\ 2 )^) 

<c'R- 2 [ |^(x)| 2 c/x. 

J{x:\xi\>R-l} 

From this (|3.8|) follows, and therefore Lemma 13.71 is proved. □ 

Next, we define an operator £r, which will harmonically extend func- 
tions defined on § 2 x M 3Ar " 3 (see (gUj) and lETTfll) below). 

We define by Yt >m {u) the normalised (in L 2 (§ 2 )) real valued spherical 
harmonics of degree £, £ G N , with m — 0,1, ... , 2£+l. Then {Yi jm }i tTn 
constitutes an orthonormal basis in L 2 (§ 2 ); they are the eigenfunctions 
for £ 2 , the Laplace-Beltrami operator on S 2 (£ 2 Y^ m = £(£ + l)Yi m ). 

Since L 2 (§ 2 x M 3A ^ 3 ) = L 2 (§ 2 ; L 2 (R 3N ~ 3 )), we have, for G L 2 (§ 2 x 
M 3Ar ~ 3 ), 

oo 2£+l 

0(^, Xl ) = J2J2 Y * 

,m(^ l -')0£,m(-^l) ) 

f=0 m=0 

0£.m(xi) = / i/ >m (w)^(a;,xi)dw. 

Note that, with T the Fourier-transform on L 2 (M. 3N ~ 3 ), 

^ m (x0 = - 1 / e^- 1 (^, m )(k 1 ) dki . (3.12) 
V2vr diV - d jR3iv-3 



/ 2^) 37V ~ 3 JR3N- 

Define 

[5 fl 0](r,w,xi) (3.13) 

oo 2£+l „ 

= TTr^FT E E / 
(V27r) 3iV - 3 ^^ 7 R 3iv-3 

where / = f e £ fl is the solution (given by Lemma IA.1I in Appendix EJ 
to the equation 

r + V-[^±^ + k 2 ]/ = 0, re^+oo), (3-14) 
satisfying /(#) = 1, \f(r)\ < 1 for all r > R. 
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Using again Lemma lA. II we find that 

5fl0eL°°([i2,oo);L 2 (S 2 x R 3N - 3 )). 

Furthermore, by the definition of f e ^ R , £ R (fi satisfies (in the sense of 
distributions) 

A£ R (f> = in {x G R 3N \ [x^ > R}. (3.15) 
We also have that 

T R E R (j> = 0. (3.16) 
By analogy with Q(Ri, R 2 ) from (|3.2j) . we define 

Q(Ri, R 2 ) = {x G R 3N \R 1 <\x 1 \ < R 2 }. (3.17) 
With this notation we have the following L 2 - bound: 
Lemma 3.9. For all R > and <p G L 2 (§ 2 x M 3Af ~ 3 ) ; 

H^^llL 2 (f2(R,3i?)) - 3 ^ 3/2 H0lU 2 (S 2 xlR 3JV - 3 )- (3.18) 

Proof. Using the definition (see ([3.13)) ) of £ R , the Plancherel theorem, 
and Lemma fA. 11 we get 

/ ( / / |£r<K?~, oj, Xi)| 2 duj dki Jr 2 dr 

= (27T) 3 - 3Ar / / Ir^Ml 2 !^^)^)! 2 

J R 

This implies the conclusion of the lemma. □ 

Lemma 3.10. Let ip be an eigenfunction of HN lt N 2 with eigenvalue E 
and assume that E < inf cr css (/fjv li jv 2 ). 

Then there exist constants Ro,c > such that for all R > Ro and 
allj e{l,...,N}, 

POO 

/ pj(r)r 2 dr ^cRZpjiR). (3.19) 
Jr 



In particular, 



POO 

/ p(r)r 2 dr <cR 3 p(R). (3.20) 
Jr 
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Proof. Clearly ()3.2()j) follows from ()3.19j) by summation over j, so we 
will only prove (|3.19|) . Without loss of generality, we only consider 
j = 1 and therefore aim to prove 

/ pi(r)r 2 dr < cR^p^R). (3.21) 
Jr 

We define (with the previously defined operators £r, T r , see (13. 7|) 
and (jniH) 

u = £rT r iI), 

as a function on R 3Ar . The function u does not necessarily satisfy 
the antisymmetry properties from Q(ifjVi,Ar 2 )- Therefore, denote, for 
x G R 3 , by m x the function on R 3JV_3 defined by 

u x (x 2 , . • • , x N ) = u(x, x 2 , ■ ■ ■ , x N ). 

We stress that u x is not a derivative of u. With this definition u x has 
the useful symmetry property 

u x e Q(H Nl _ 1>N2 ) . (3.22) 

From Lemma f3. 91 we get the inequality 

/ |n(x)| 2 rfx < Q^llT^ll^^x^-a) = 9R 3 p 1 (R). (3.23) 

JQ{R,3R) 

Define H Nl)N2 (R) as the operator obtained by restricting H to the space 
7~Ln 1 ,n 2 {R) '■ = 

w 2 > 2 (r 3 \WJV) n w^ 2 {r 3 \ WJ*j)) ® ^f-i,iv 2 (K 3Ar " 3 ) • 

That is, we impose Dirichlet conditions at radius R on the first electron 
coordinate, and the symmetry conditions on the last N — 1 electron 
coordinates. 

Let tp e H Nl , N2 (R) be normalised (in L 2 (R 37V )). Then, since \xi\ > 
R, it follows that 

N z 

(<p,H NuN2 (R)<p) >{i P ,[Y i (-A j - 



i=2 



7?" 



By the HVZ-theorem (see jlll Theorem XIII. 17']) the term in ( ■ , ■ ) 
(on the right side) is bounded below by inf a^ B (H Nl)N2 ). Thus, for any 
e > there exists R' > such that for all R> R', 

inf cr(H NliN2 (R)) > inf a ess (H NljN2 ) - e. 



POSITIVITY AND DECAY OF THE ELECTRON DENSITY 17 

Since, by assumption, E < inf 0"e S s(-ffjVi,Jv 2 )) the operator Hn u n 2 (R)—E 
is invertible for R sufficiently large, i.e., for all R> Rq for some Rq > 0. 
Let C e C°°(R), < C < 1, be a function such that 

C(t) = 1 for \t\ < 2, £(t) = for \t\ > 3. (3.24) 

With ( as above and R > 0, we define (r{xi, . . . ,xn) '■= C{\ x i\/R)- Let 
v := ip—( R u. Then T R v = 0, so we see using (|3.22|) that v G T~tNi,N 2 {R)- 
A calculation gives 

(-A + V - E)v = -(V - E)( R u + 2VCrVm + (A( R )u. (3.25) 

Since H]sr lt N 2 (R) — E is invertible and v G T-Ln 1 ,n 2 {R)i we find 

v = (H NuN2 (R) - E)-^- (V - E)( R u + 2VCrVu + (A( R )u). 

It is easy to see, using that V is relatively bounded with respect to the 
Laplacian, ()3.23|) . and the support properties of ( R , that there exist 
c, d such that 

INL^x:^^!}) - HPllia^xi^lajKaii}) — c ^ Pi{R)- (3.26) 
Combining (|3.23|) and (|3.26|) we get 

IMIz3({x:.R<|:ci|}) = IICR' u + ' u llL2({x:fl<|a! 1 |}) 

- 2 (llCi?w|lL2({ x . /?< | :ri |}) + ll^ll^^x.^i^i})) 

< cRZptiR). 

This is the inequality (|3.21j) . The proof of Lemma 13.101 is therefore 
finished. □ 

The estimate ()1.12|) follows from ()3.9|) and ()3.5|) . The lower bound 
()1.13|) clearly follows from Lemma f3 . 1 01 upon inserting ()3.6|) . This fin- 
ishes the proof of Theorem 11.41 

Appendix A. 

Lemma A.l. For all £ e NU{0} ; all R > and all k > 0, the equation 

f"+ 2 -f'-[^ }1 + - 2 }f = 0, f(R) = l (A.l) 

has a solution f vanishing at infinity and satisfying 

< 1 for allr > R. 

Proof. Actually, if / is a solution of (|3.14|) . then rf is a Whittaker 
function (see [T] for details). This implies that a solution / exists 
vanishing at infinity. 
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Define, for x G M 3 , u(x) := /(|x|), then u satisfies 

-Au + Wu = 0, A\ x \= R = l 

with W(x) = e -^± + k 2 > 0. By Kato's inequality [TU1 Theorem X.27] 

we get 

-A\u\ + W\u\ < 0. (A.2) 
Let v K , k > 0, be the function on M 3 \ {0}, 

, , R_ 

\x\ 



v K (x) = —e-W-R). 



Then 

— Av K + k 2 v k = 0, v K \, ,_ R = 1, and vJx) < 1 for \x\ > R. 



\\x\=R 



So 



(-A + W)(v K -\u\)> £ -^^-v K >0, 

\x\ z 

K-N)| N=Jl = o. 

The maximum principle (see e. g. [5, Theorem 8.1]) implies (since / 
and v K vanish at infinity) that for all |x| > R, 

\u(x)\<v K (x). (A.3) 

This implies the statement of Lemma IA.1I □ 
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